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Let £ be a planet modeled by a homogeneous rigid ellipsoid with symmetry axis NS, with polar inertia 
moment J 3 , equatorial moment J and mechanical flattening coefficient r\ = (J 3 — J)/J 3 . 1 

The planet center of mass T is supposed to revolve on a keplerian orbit t — > fr(i) about a focus S: the 
orbit plane will be called the ecliptic plane and k will denote its unit normal vector {celestial north) which 
sees the planet revolving counterclockwise. 

The longitude At of Tt on the ecliptic will be reckoned from the aphelion of the ellipse; hence Ax = is 
the aphelion position, i.e. where Tt = \r*x\ is maximal: rx(0) = a(l + e), a being the major semiaxis of the 
keplerian ellipse and e its eccentricity. 

With the above conventions, rr and Ay are related by: = \rr\ — p- (1 — e cos At) -1 , p = a(l — e 2 ) and the 
Kepler's laws imply that if A is the keplerian average anomaly then: A = At + 2e sin At + |e 2 sin 2At + . . .. 

and the motion is A — > A + u>Tt, where 2it/ujt = 27ra 3 / 2 <; JV 1 ^ 2 is the year of the planet, gjy = k(ms + Tot) if 
k is Newton's constant and rar, tos are the masses of the planet and of its star. 
The planet is described by means of the following coordinates: 

1) The total angular momentum, or spin, M and its projections K = Mcosz on the axis orthogonal to the 
ecliptic ("celestial north-south axis"). 

2) The angles 7, <p, A where 7 is the angle between an axis fixed on the ecliptic (Aries line) and the spin- 
ecliptic node (i.e. the intersection of the ecliptic and the plane orthogonal to the spin: equinox line); <p is the 
angle between the spin ecliptic node and the analogous spin-equator node; and A is the average anomaly of 
the planet, which rotates uniformly at angular velocity ujt- 

It is a well known theorem in classical mechanics (Andoyer-Deprit) that the pairs (M, if) and (K, 7) are 
canonically conjugate variables (see [G], p. 318) and, introducing an auxiliary variable B canonically conjugate 
to A, the energy of the system can be written: 



One could avoid introducing B: but one would then have a time dependent hamiltonian, which I do not like. 

The complete description of the rigid body configuration would require an extra pair of canonically conju- 
gated varibles, namely (L, ip) where ip is the angle between a comoving axis fixed on the equator and the 
spin-equator node and L is the projection L = M cost? of the spin on the NS axis. However, by symmetry, 
the hamiltonian does not depend on ip; hence L is a constant of motion and one can see that it only gives 
an additive contribution to the energy, which is dropped in (1). 

I shall consider the following approximation for V: 

a) the integral is developed in powers of the ratio R/a and the expansion will be truncated to the order 2 
included neglecting the orders 4 and higher. 2 

b) the expression thus obtained will be developed in powers of the eccentricity and the expansion will be 
truncated to 2d order. 

The model thus obtained will be called the D' Alembert precession-nutation model. The e = model was 
in fact used by D'Alembert, [L], to deduce his celebrated theory of the equinox precession for the Earth and 
the theory of the lunar plane precession. 

It is not interesting to limit ourselves to the D 'Alembert 's case e = because the phenomenon of large 
variations of the inclination axis are in this approximation not possible. 

Define ujd = M / J , calling it conventionally the daily rotation, and consider the line C in the plane of the 
action variables A = (M, K) defined by: 
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1 this means that if _R is the equatorial radius then the polar radius is R/(l + 2rj) 1 / 2 . 

2 n 

only the 2d order is non trivial among the first three orders. 




(1) 



M = Jlujj — 2Jujt = const, 



K < K < K x 



(2) 
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and call ^4° and A 1 its extremes. Since K = Mcosi one realizes that proceeding on the above line from 
A to A 1 the spin inclination angle changes between i and i\ and the system is locked in a 1 : 2 spin orbit 
resonance. 

A motion in which the projection in the (M, K) plane closely follows the line £ is therefore a motion during 
which a variation of the inclination axis is observed. Along such line, which will be fixed in the following, 
the ratio cosd = L/M is a constant which will therefore also be fixed: and which will be supposed different 
from ±1, (i.e. the angle d, nutation constant, between the spin and the NS axis is supposed different from 
or 7r). 

The eccentricity will be taken e = rf for some c > 1 so that setting r\ small implies that the eccentricity is 
also (much) smaller. 

The question that will be investigated, see [CG], is whether there exists a trajectory X e ^(t) starting at 
t = in j4e,r; — {M e ^ , -Ke,?j) near A , with phases (7 6j7; , <^ e;J? , A ej7? ), and reaching after a long enough time 
the vicinity of A 1 . 

If such a trajectory exists for all n small enough, but different from 0, one says that there is a Arnold 
diffusion phenomenon. 3 

What is remarkable, when the diffusion phenomenon happens, is that it happens in spite of the fact that 
K, M are adiabatic invariants, i.e. they are variables which in any averaging approximation could only change 
by quantities which are infinitesimal as e, n — > 0: while in the above situation K or i change by a quantity 
~ i\ — io independent on e,r]. 

The existence of Arnold diffusion says therefore that no averaging approximation can be even approximately 
correct over very long time scales or for infinite time. 

This is not at all a knell for the widely used averaging methods: but it shows the interest of estimating the 
time scales over which the phenomenon happens: it is becoming well understood that the correctness of the 
averaging methods give very accurate results over time scales that are often of the order of the age of the 
universe, [CeG]. 

The result, [CG], on the above question, in the case of the described precession nutation model is the 
following. 

Theorem: Suppose that [io,ii] does not contain the values 0, -|,7r and suppose that the angle t? between the 
NS axis and the spin axis does not vanish in the configurations in £ line. There are constants c, b, d such 
that if e = rf then for all r] small enough, but non zero, there are phases = (7^, <f n , X v ) and actions 
A n = (M rj ,K ri ) such that: 

1) the trajectory X v (t) = (A(t), $(t)), starting in the configuration A(0) = A v , $(0) = ^ v stays close, in the 
A-variables, within n b to the line £ during a suitable time > T v . 

2) the time T n is such that A(T V ) is close (i.e. within n h ) to A 1 while A(0) is close to A a , (i.e. within n h ). 

Therefore in the above precession problem Arnold's diffusion is possible. The angle of spin inclination can 
change by an order of magnitude 0(1) in a long enough time, no matter how small is the spin orbit coupling 
(i.e. n), provided it is non zero and provided the eccentricity is small enough compared to the flattening. 
And it becomes interesting to get an idea of the size of the time T v . 

In [CG] an explicit estimate for T v is derived: 

Theorem: in the situation of the previous theorem there is a constant f such that for n small enough the 
diffusion time T v can be estimated to be not longer than: 

T„=u,~V e V /2 r 2 ^)- 2 / «,(£)= sin min ( 1 ~ C ° si cos 2 i sin i) (3) 

n T yv ; i <i<n VI + cosi / v ; 

Hence the estimate is a very very long time, compared to the times over which one expects the averaging 
methods to have some validity, which have scale of order ~ exp frj^ 1 / 2 . Also the constants b, c, d, f, whose 
values can actually be computed, see [CG], turn out to be very poor for any practical application. 

But to get to practical applications a lot of estimates would have to be refined and it is not completely 
obvious that this is really impossible. One can recall, on this respect, that for a long time it was stated by 
some people that the KAM estimates were too bad for any practical use: and this turned out to be grossly 
incorrect, [CeG],[CC],[LR]. In the case of diffusion the situation looks much harder, perhaps desperately so: 
but we shall see. 

this refers to one possible definition of Arnold's diffusion, [A]: it is a special case of a general definition proposed in [CG]. But 
there are other definitions: different and, often, not very precise, 
it has been remarked above that the angle $ is constant on the line C. 
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This shows that the estimate of the time scale for the diffusion, i.e. essentially the time one has to wait 
to see a violation of the predictions of the averaging approximations, diverges if one tries to bridge gaps in 
inclination containing i = 0, ^, tt. 

In particular the above estimate diverges if one tries to find a trajectory in which the spin sign is reversed 
(i.e. i goes through ^). The difficulty in constructing such trajectories might be a manifestation of a physical 
phenomenon and not just a defect of the techniques of proof; it seems, indeed, that the chaotic motions of 
the planetary axes are unable, in absence of dissipative effects, to change the spin sign, i.e. the sign of K, 
[LRo] . 

Formula (3) is suggestive as it leads to the conjecture that (in general) diffusion along a path C in action 
space might really take place over a time scale dependinng exponentially on some power of the coupling 
constant times a coefficient determined my maximising a suitable functional defined on C: therefore I call 
(3), by analogy, a large deviation formula. 

I conclude with a technical comment: the analysis is based on the fact that the hamiltonian system under 
consideration has three very different time scales: namely the daily time scale (u;^, , coinciding with the year 
time scale because the free system is in a 1 : 2 resonance), the resonance time scale of order O^ -1 / 2 ^ 1 ) 
(describing the characteristic time of oscillation transversal to the resonance) and the equinox precession scale 
of order C^r; -1 ^ 1 ) (arising from an application of the method of averaging and describing the precession 
of the equinox line, when e = 0). The three time scales become very different from each other when 77 — > 0: 
this is due to the degeneracy inherent in the problem where the unperturbed hamiltonian does not depend 
on the variable K. And it implies the phenomenon of large angles at the homoclinic intersections, see [CG], 
which in turn implies the existence of diffusion. Deeply different is what happens near resonances of non 
degenerate systems, where all the action variables appear non linearly in the unperturbed hamiltonian: in 
such cases there are only two distinct time scales describing the motions near the resonances and the theory 
of diffusion is much harder (because the homoclinic angles are exponentially small, see [CG], [G2]). 

Thus the above theory is only one more instance of the (well known) fact that the degeneration intrinsically 
present in all celestial mechanics problems is very often, in fact, making the theory easier rather than harder 
(as sometimes naively claimed). 

For connections between the techniques used in the theory of the above problem and the KAM theorem see 
also [G2],[G3]. 
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